A New Regime of Anomalous Penetration of Relativistically Strong 
Laser Radiation into an Overdense Plasma 

A. Kim, 1 F. Cattani, 2 D. Anderson, 2 and M. Lisak 2 
institute of Applied Physics, Russian Academy of Sciences, 603600 Nizhny Novgorod, Russia 
2 Department of Electromagnetics, Chalmers University of Technology, S-412 96 Goteborg, 

Sweden 

Abstract 

It is shown that penetration of relativistically intense laser light into an overdense 
plasma, accessible by self-induced transparency, occurs over a finite length only. The pen- 
etration length depends crucially on the overdense plasma parameter and increases with 
increasing incident intensity after exceeding the threshold for self-induced transparency. 
Exact analytical solutions describing the plasma-field distributions are presented. 
PACS number(s): 52.40.Nk, 52.35.Mw, 52.60. +h, 52.58.Ns 

In the past few years there has been much research devoted to the nonlinear interaction of 
superintense laser pulses with plasmas At intensities where electrons quiver with relativistic 
velocities, the interaction can be characterized as nonlinear optics in relativistic plasmas, and 
new regimes, not evident at nonrelativistic intensities, may appear. As was previously shown, 
superintense electromagnetic radiation can propagate through a classically overdense plasma 
due to the relativistic correction to the electron mass, the so called induced transparency effect, 



|2|]-|[TC|. The present work has resulted in the identification of a new fundamental process in 
relativistic laser overdense plasma interaction. 

In order to understand the nonlinear regime of interaction of superintense laser light with 
an overdense plasma, it is enough, without loss of generality, to consider a stationary model. 
We present here a new class of exact analytical solutions describing the penetration of an 
electromagnetic wave normally incident onto a cold, overdense plasma with a sharp boundary. 
In particular, we show that, when the incident intensity exceeds the threshold for self- induced 
transparency, the laser energy penetrates into the dense plasma without any losses, but over a 
finite length only. At the same time, the electron density distribution becomes structured as a 
sequence of electron layers separated by about half a wavelength wide depleted regions, so that 
this strongly nonlinear plasma structure acts as a distributed Bragg reflector. 

The ultrahigh intensity laser-plasma interaction is described by the relativistic equation 
of motion and the equation of continuity for the electrons together with Maxwell's equations. 
Ions are treated as a uniform neutralizing background. We will consider circularly polarized 
laser radiation with normalized amplitude of the vector potential eA/mc 2 = (a/y/2)Re[(y + 
iz)exp(iut)] normally incident from vacuum (x < 0) onto a semi-infinite plasma (x > 0). 
Assuming a stationary regime, the basic equations may be written in the form 

0" = n o (n-l), (1) 

a" + (1 - — n)a = 0, (2) 

7 

(p' = 7' if and only if n(x) ^ 0. (3) 

(where variables are normalized as x =>■ ujx/c, n =>- n/n , n Q = uj 2 /uj 2 , uj is the carrier frequency 
of the laser radiation, u p is the plasma frequency of the initial unperturbed plasma, 7 = 
(1 + 2 ) 1 / 2 is the relativistic factor, n and <p are normalized electron density and scalar potential 
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of the plasma respectively). Eq. (3) indicates that only in the region where the electron density 
n(x) 7^ the ponderomotive force, j'(x), must be compensated by the force of a longitudinal 
field. This statement will be important for constructing solutions of interest in the present 
analysis. 

For homogeneous ion density, the system has the following Hamiltonian 
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He = 2(l + a*) a ' 2 ~ 2 {2noVYT ^ ~ ° 2) (4) 

which was analyzed in ||. 

As we are interested in a semi- infinite plasma, we first consider the case when n(x) — > 1 
and both a(x) and a'(x) vanish as x —>■ oo. In this case the integral of motion equals 

He = -n = Ho- (5) 
Eq. (H) can then easily be integrated to yield the following single-parametric solitary solution 
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a ( X ) = ,,2n. ../nun ! i ( 6 ) 



where e$ = 1 — n is the dielectric permittivity of the plasma and the parameter x^ defines 
the position of the maximum of the function (0) which is given by A m = 2[n (n — l)] 1 ^ 2 . 

For n a > 1.5 this solution contains a region where the electron density is negative, which is 
clearly unphysical. Requiring the minimum electron density n m = 1 — 4(n G — l) 2 to be positive, 
i.e. n m > 0, we obtain a condition on the background plasma density and an argument in 
favor of the depletion region: for n > 1.5 we have to take into account only that part of the 
solution where the corresponding electron density is positive, n(x) > O.By using the part of the 
solution (H) corresponding to a ponderomotive force that pushes electrons into the plasma and 
matching it to the vacuum one, an exact expression for the intensity threshold of self-induced 



transparency was found, |TTJ. It should be noted that the Hamiltonian (gj) corresponds to a 
zero energy flux, but there are solutions with non-zero flux as well, |2|, |J . They can arise even 
in correspondence to incident amplitudes smaller than the penetration threshold. However, for 



the realization of such nonlinear hysteresis-like solutions, ]T0| , |T2] , preliminary modifications of 
the plasma must be induced by extremely intense fields. For the description of the steady-state 
of a problem involving laser radiation with turn-on shape it is natural to choose the zero flux 
approach. 



The left side of the function gives rise to an unbalanced ponderomotive force pulling 
electrons out of the plasma towards the incident wave and, at first glance, it seems that the 
charge quasineutrality condition cannot be satisfied. However, after a length of about half a 
wavelength, these electrons will be stopped by the ponderomotive force acting in the opposite 
direction. Thus, in the general case, we may expect that the full plasma-field structure will 
consist of a sequence of alternating depletion and nondepletion regions. This can be understood 
from Fig. 1, where the phase portraits described by Eqs (l)-(3)and by the corresponding 
equation for the depletion (vacuum) region are presented. For a half-space plasma the limiting 
case corresponds to a motion with infinite period along the separatrix determined by Eq. fl5|), 
i.e., an exponentially decreasing field inside the overdense plasma. Going backwards towards 
the initial vacuum-plasma boundary, before the last semi-infinite electron layer there must exist 



2 



a depletion region. Here the amplitude of the field corresponds to a forward going wave along 
the incident direction with an intensity below the threshold value (corresponding to the motion 
along the first circular trajectory in phase space, coming from the vacuum Hamiltonian) . Then, 
in front of this layer we have to put another electron layer, where the solution for the field a(x) 
follows from the Hamiltonian (Q) with a magnitude 7i > —n (corresponding to an oscillatory 
motion about zero). This construction is repeated until the initial plasma boundary is reached. 
At the boundary between each depletion and nondepletion region, the solution for the field 
must satisfy continuity conditions for both a and its derivative a'. It is also clear that there 
exists a family of stationary solutions that differ from each other by the number of electron 
density layers and their shapes. When the incident amplitude is increased, the number of layers 
will increase as well, as follows from the phase portraits in Fig. 1. 

In order to quantify the above discussion, we present a more rigorous analytical description. 
Starting from inside the plasma region, the solution for the field is an exponentially decreasing 
function of the spatial coordinate, in fact a part of the localized solution given by Eq. (||) with 
He = 7~Lo for x < x < oo. The point x = x can be determined self-consistently from the 
global solution and the boundary conditions 

a(x = x ) = a , a'(x = x ) = a' (7) 

The next region (x\ < x < xq) must be a depletion layer where n(x) = 0. The second 
boundary position of this depletion layer, x±, must also be determined self-consistently. The 
Hamiltonian here is the vacuum one and reads 

n v = ~ (a' 2 + a 2 ) = ~ K 2 + a 2 ) = H x (8) 

where we have taken into account the boundary conditions for the field and its first derivative. 
The solution of the field in this depletion region reads 

a(x) = Ai cos(x + if)x). (9) 

Using ([?D we have: A x = (a(, 2 + a 2 ,) 1 / 2 = (2H1) 1 / 2 , ip\ = — arctan(do/a ) — ^o- The boundary 
position Xican be calculated by integrating Poisson's equation within the interval X\ < x < x : 

<f>'(x ) - (f}'(xi) = -n o (x - xi). (10) 

which, together with relation (3), leads to a transcendental equation for £ = Xq — X\\ 

£ = 9(0-9(0) (11) 

where 

g(0= 2^rTA!^TM^ ' 6 = arcta H^j ( ' 

Since the left-hand side of Eq. flTT| ) is a linear function while its right-hand side is a periodic 
function of £, this equation has a non-trivial solution only if <?'(£ = — £1) > 1, i.e. 

However, if n Q < 1.5, this condition is never satisfied, because its maximum value reaches one at 
n Q = 1.5. Consequently, for plasma densities n Q < 1.5 we conclude that there are no stationary 
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regimes of anomalous penetration: there can only exist dynamical solutions. For n > 1.5, 
Eq. (pTT|) always has non-trivial solutions, which can be found numerically. Formally, Eq. (|TT| ) 
admits several roots but we have to consider only those that correspond to a positive electron 
density. 

Having solved Eq . (|TTD we know a\ = a{x\) and a[ = a!{x\). The next region must again be 
an electron layer. Denoting this region as x 2 < x < x\, the Hamiltonian (|j) will here be 

n E = n E (a 1 ,a' 1 ) = n 2 . (14) 

with the corresponding field solution written in terms of two-parameter elliptic functions for 
the field strength as 



a[x) 



2q C n{ey\x- x^))\ _ 

Ii727 n 2 < n o, 



2+[( 9 2 +l) 1 /2_l] sn 2[ £ V- ! ( :E _ 2 .(2))]> * V 
2gsn[((£2+l) 2 -ng) 1 / 2 (a:-x( 2 ))/2] ' 
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2 -sn2[((e 2 +l)2_„2)i/2( ;c „ ;! ,(2))/2] 



7^2 > "-o 



The two parameters, £2 and are determined by the boundary conditions at x — x\. Here 
e 2 = (n 2 + 1 + 2ft 2 ) 1 /2, q = [( £2 + n ) 2 - l] 1 ^ and g = [( £a + n + l)/(e 2 + n - l)] 1 /*, 
fc = [( n 2 - (e 2 - 1) 2 )/4£ 2 ]V2 an d k = [(e 2 - l) 2 - n 2 )/(e 2 + l) 2 - n 2 )] 1 ^ are the moduli of the 
elliptic integrals of the first kind respectively for the two cases. This solution is realized over 
about a quarter period, where the electron density is positive and is given by 



n ( x ) = Vl + ( H 2 + -n VT+^ - a 2 ) (16) 
n Q \ 2 / 

The second boundary position of this layer can now be defined with a certain arbitrariness, 
namely, it can be taken within the whole interval x* < x 2 < x** , where x* is such that n(x*) = 0, 
while x** is defined by the existence of a solution of the transcendental equation (0). This 
problem comes up for every layer, therefore the global solution, for a given incident amplitude, 
in general may not be unique because the thickness of each plasma layer is not uniquely fixed, 
the realization of each specific solution depending on the prehistory of the process. This is also 
confirmed by our preliminary simulations based on a hydrodynamics approach where, at fixed 
incident amplitude, the thickness of the various electron layers in the quasi-stationary stage is 
quite different depending on the turn-on shape of the incident radiation. The appearance of 
non- unique solutions is typical of phenomena involving nonlinear media (see, for instance, ||12|Q . 

For the sake of concreteness, in what follows, we choose the next boundary position x 2 as 
the point where the electron density vanishes, i.e. n(x 2 ) = 0. In this case, the transcendental 
equation ( [11]) always has nontrivial solutions because, if the condition given by Eq. ( |T3| ) 
is satisfied for the semi-infinite layer, it will automatically be satisfied for all the precedent 
electron layers, for the point where the electron density vanishes. 
As follows from (|16D, the field at this boundary is such that 

3n (l + a 2 2 y/ 2 = 2(a 2 2 -H 2 ) (17) 

and, making use of the equation for the Hamiltonian (|T4"D, we can calculate the first derivative of 
the field. Thus, for the next depletion layer located within £3 < x < x 2 , we know the boundary 
values of a 2 and a' 2 that are required to generate a solution by using the same procedure as 
before for x\ < x < x . This procedure is repeated for every layer, x i+ i < x < x i7 up to the 
final layer and if Tii > n , the field solution has to be taken as in the second expression of flT5|). 
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The last layer will be a depletion region < x < xn where x = is the real plasma- vacuum 
boundary. At this point the electric field due to charge separation must vanish. Integrating 
Poisson's equation over the plasma interval, the plasma neutrality condition gives 



which defines the last free parameter. Thus, matching the solution for the field to the vacuum 
solution at the electron layer boundary x = xn we can relate the obtained plasma- field distri- 
bution to the incident electromagnetic wave. In this way we can construct an exact stationary 
solution for the anomalous penetration regime. An example of such a solution is presented in 
Fig. 2. Notice that the maximum of the electron density in a layer increases from layer to 
layer, reaching an absolute maximum in the last layer nearest to the vacuum boundary while 
the width of the layers become more and more narrow. This may easily be understood from 
Fig. 1; higher lying trajectories have higher values of the integrals of motion. It also means 
that at higher incident wave amplitudes there are more layers, so that the penetration length 
will increase with increasing incident intensity. These solutions can naturally be extended to 
cover the case of plasma slabs with finite thickness, if the penetration length is smaller than 
the slab thickness. 

The following consequence of the previous analysis should be emphasized: For a fixed am- 
plitude the laser field penetrates only a finite length into the overdense plasma. Consequently, 
if the incident laser pulse has a finite duration, the electromagnetic energy deposited in the 
plasma will be reflected back into the vacuum after the laser field has vanished. Evidently the 
transient regime will be more complicated as vacuum (depletion) regions surrounded by elec- 
tron layers will show a resonator-like behaviour, with the electromagnetic energy being excited 
by the incident pulse. We expect this field structure to evolve and to be slowly reflected back 
into the vacuum region. 

In conclusion, we have presented a new class of exact analytical stationary solutions describ- 
ing a new feature of the interaction between a super-intense laser and an overdense plasma. 
This analysis shows how, depending on the initial plasma density, the interaction might result 
in the generation of a new plasma-field structure consisting of alternating electron and "vac- 
uum" regions, with the electromagnetic energy penetrating into the overdense plasma over a 
finite length determined by the incident intensity. 

This work was supported in part by the Russian Fund for Fundamental Research (grants 
No. 98-02-17015 and No. 98-02-17013). One of the authors (F.C.) acknowledges support from 
the European Community under the contract ERBFMBICT972428. 

Captions 

Fig. 1 Phase portrait of the system for n a = 1.6 and homogeneous ion density, following 
from Eqs. (|J) and (|8|), see also Fig. 2. The dashed lines stand for vacuum regions, the con- 
tinuous lines stand for plasma regions, the actual trajectory is given by the thick line and it 
runs clockwise. The large dashed line denotes the regions where the electron density is negative. 



Fig. 2 The continuous line represents the plasma-field structures in a semi-infinite plasma 
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initially occupying the region x > 0, for n a = 1.6, the unperturbed electron density being 
represented by the dashed line. The dotted line represents the resulting electron density distri- 
bution. All quantities are dimensionless. 
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Fig. 2 




